Draft version February 2, 2008 

Preprint typeset using I£T^X style cmulatcapj v. 7/15/03 



^1- 
O 
O 

(N 

Oh- 
< 



> 

o 
o 

Of 

6 

c3 



X 



LIMITS ON /-BAND MICRO VARIABILITY OF THE GALACTIC BULGE MIRAS 
P. R. Wozniak, K. E. McGowan, AND W. T. Vestrand 

Los Alamos National Laboratory, MS-D436, Los Alamos, NM 87545 
Draft version February 2, 2008 

ABSTRACT 

We search for microvariability in a sample of 485 Mira variables with high quality /-band light curves 
from the second generation Optical Gravitational Lensing Experiment (OGLE-II). Rapid variations 
with amplitudes in the ~0. 2-1.1 mag range lasting hours to days were discovered in Hipparcos data 
by de Laverny et al. (1998). Our search is primarily sensitive to events with time-scales of ~1 day, 
but retains a few percent efficiency (per object) for detecting unresolved microvariability events as 
short as 2 hours. We do not detect any candidate events. Assuming that the distribution of the event 
time profiles is identical to that from the Hipparcos light curves we derive the 95% confidence level 
upper limit of 0.038 yr -1 star -1 for the rate of such events (1 per 26 years per average object of the 
ensemble). The high event rates of the order of ~1 yr -1 star -1 implied by the Hipparcos study in the 
Hp band are excluded with high confidence by the OGLE-II data in the / band. Our non-detection 
could still be explained by much redder spectral response of the / filter compared to the Hp band or 
by population differences between the bulge and the solar neighborhood. In any case, the OGLE-II 
/-band data provide the first limit on the rate of the postulated microvariability events in Mira stars 
and offer new quantitative constraints on their properties. Similar limits are obtained for other pulse 
shapes and a range of the assumed time-scales and size- frequency distributions. 
Subject headings: stars: variables, AGB and post-AGB, activity 



1. INTRODUCTION 

Mira variables are a subclass of Red Variables, also 
known as Long Period Variables (LPVs). They are low 
mass AGB giants, mostly with M type spectra (e.g. 
Gautschy & Saio 1996). Their large amplitude light vari- 
ations on time scales of ~1 year are thought to arise due 
to the pulsation instability (e.g. Keeley 1970). There are 
intriguing reports of rapid light variations in Mira stars 
on time scales of a few days (Maffei & Tosti 1995) or 
even hours and minutes (Schaefer 1991) with amplitudes 
in the range 0.1-1.4 mag. The most systematic study 
of microvariability in Mira stars is based on Hipparcos 
data (de Laverny et al. 1998, hereafter dL98). It reports 
the discovery of 51 events in 39 out of 239 observed Mira 
stars. The amplitudes of events detected by Hipparcos 
are between 0.2 and 1.1 mag. To best of our knowledge, 
so far there are no published estimates for the rate of 
such events. Other reports of rapid variations in broad 
band photometry of specific Mira variables include Smak 
& Wing (1979) or more recently de Laverny et al. (1997) 
and Teets & Henson (2003). From existing data very lit- 
tle can be inferred about the spectral properties of those 
rapid variations, but microvariability in several spectral 
features has been reported (Odell et al. 1970, Kovar et 
al. 1972). On the other hand, Smith et al. (2002) do not 
detect any rapid variations in the COBE DIRBE mid- 
infrared light curves of 38 Miras. 

Willson & Struck (2002) explored possible interpreta- 
tions and favored one involving hot flashes due to inter- 
action of the extended Mira atmosphere with a Jovian 
planet or a brown dwarf. This process was studied nu- 
merically by Struck, Cohanim & Willson (2002, 2004). 
The scenario can explain the time scale of the observed 
microvariability and produces energies and luminosities 



roughly of the right order of magnitude, but it ignores the 
fact that dL98 have seen approximately equal numbers 
of positive and negative flux excursions. 

Here we use a large set of high quality /-band photom- 
etry from the Optical Gravitational Lensing Experiment 
(OGLE) to constrain the rate of microvariability for Mira 
stars found in the Galactic Bulge. 

2. DATA AND SAMPLE SELECTION 

2.1. OGLE-II light curves 

We use the photometric data collected, processed and 
published by the second generation OGLE experiment 
(OGLE-II; Udalski, Kubiak, & Szymanski 1997). Our 
analysis is based on photometry from Difference Image 
Analysis (DIA) of the images collected during observ- 
ing seasons 1997-1999 (Wozniak et al. 2002, Wozniak 
2000). The DIA technique relieves numerous complica- 
tions associated with crowded fields and enabled routine 
photometry as accurate as 0.5% for bright unsaturated 
stars like the brightest AGB variables in the Galactic 
Bulge. This data set contains /-band light curves for ap- 
proximately 2xl0 5 candidate variable objects with the 
average time sampling of once every 1-3 days and typi- 
cally 250 measurements per object. We defer a detailed 
discussion of the time sampling to Section l3~Tl Observa- 
tions span the magnitude range between 10.5 and 19.5. 
The total time baseline of the data set is 3 years with 2 
inter-season gaps, each lasting about 4 months. 

2.2. Selecting Mira variables 

The General Catalog of Variable Stars (GCVS; 
Kholopov 1998) states that Mira variables have F-band 
light amplitudes between 2.5 and 11 mag, well pro- 
nounced periodicity and periods ranging from 80 to 1000 
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Fig. 1. — Period amplitude distribution (left) and histogram of 
amplitudes (right) for a sample of red variables from OGLE-II data. 
The location of the vertical line indicates the adopted amplitude 
threshold for Mira variables. 

trum. The histogram of amplitudes for LPVs with log(P) 
larger than 2.0 is bimodal (Mattei et al. 1997) with rea- 
sonably good separation between lower amplitude semi- 
regular variables (SR) and Miras (M). 

We begin by calculating the total amplitude A for all 
2xl0 5 objects in the OGLE-II database of DIA photom- 
etry using a running median filter with 11 points and 
taking the difference between extreme values. A sample 
dominated by LPVs is obtained by requiring A > 0.8 
mag. This selects 842 stars variable on time scales of 2 
weeks or longer. Periods are obtained using the Analy- 
sis of Variance algorithm (AOV; Schwarzenberg-Czerny 
1989). The period-amplitude diagram and the histogram 
of amplitudes for all 842 stars is shown in Figure ^ The 
gap separating semi-regulars and Miras is clearly visi- 
ble at A ~ 1.2 mag (c.f. Cioni et al. 2001 and Eggen 
1975). As an independent cross-check, we examine near- 
infrared colors of candidate stars from the Two Micron 
All Sky Survey (2MASS 1 ). We require a positive identifi- 
cation with a 2MASS point source object within a 1" ra- 
dius of the OGLE position. Almost all identified stars 
have (J — H) and (H — K s ) colors consistent with AGB 
variables observed through the reddening screen equiv- 
alent to A v ~ 0-2 mag of visual extinction. We reject 
only 3 out of 705 positively identified stars, those having 
(J -H) < 0.75 and (H - K s ) < 0.3. After verifying the 
2MASS colors we select a sample of Mira variables using 
amplitude (A > 1.2 mag) and period (P > 100 days). 

2.3. Quality thresholding 

Our analysis is limited to the highest quality data, i.e. 
light curves with the faintest measurement of / < 15.0. 
In this magnitude range the photon noise is not a domi- 
nant part of the error budget and the influence of crowd- 
ing is relatively small. Further, only measurements with- 
out processing flags are accepted to avoid possible prob- 
lems due to saturation, unreliable PSF fits, exceptionally 
large seeing FWHM etc. Finally, measurements with un- 
usually large error bars are rejected. The median and the 
standard deviation a around the median are calculated 
for all error bars in a given light curve. Points with an 
error bar more than 4<r larger than the median error are 
removed from the analysis. The latter cut removes only 
465 points out of 105,890 and has no influence on our 
conclusions ( Section I3.1fl . The final sample consists of 
105,425 /-band measurements of 485 Mira type variables. 

3. METHOD 

This publication makes use of data products from the Two 
Micron All Sky Survey, which is a joint project of the Univer- 
sity of Massachusetts and the Infrared Processing and Analysis 
Center/California Institute of Technology, funded by the National 
Aeronautics and Space Administration and the National Science 
Foundation. 



3.1. Signature of microvariability 

If microvariability events by tenths of a magnitude or 
more with durations larger than ~0.1 days are present, 
then they should be visible as a tail of photometric out- 
liers. Some of the outliers could be statistical fluctuations 
in the measurements, so by assuming that all such out- 
liers are due to rapid variability in the sources we obtain 
an upper limit on the frequency of such events. One can 
only afford such a conservative assumption and still de- 
rive a useful limit when the data quality is very high and 
the contribution from random measurement errors to the 
tail is negligible. 

To define the normal behavior of an object with respect 
to which we measure the outliers, we fit a smooth curve 
to the data. Using x 2 minimization, each season for the 
light curves is fitted separately with a model consisting 
of a set of harmonics: 

7 

mi(t) = A + (Ak sin(fcwt) + B k cos(fcwt)), (1) 

where u> = 2ir yr _1 . This model is smooth on time scales 
of 26 days or longer. The choice of any particular model 
is not important as long as it is smooth on time scales 
of about 1 month and still fits the long term behavior of 
the Mira flux. Oscillating functions can usually achieve 
that with fewer parameters than other basis functions. 

It is important to verify that the model does not over- 
fit the data, as this would suppress the detection of out- 
liers. In Figure [3 we present a few examples of OGLE- 
II light curves for Mira stars used in our calculations. 
There is some limited potential for over-fitting near the 
beginning and the end of each observing season, espe- 
cially if one or two points are somewhat separated from 
the rest of the measurements in a particular season. This 
is unavoidable as the model for normal behavior of the 
flux near such loose points is effectively undefined. But 
within the observing season there is very little possibil- 
ity that a model smooth on time-scales of 1 month (Fig- 
ure can fit individual points. The temporal spacing of 
measurements is too fine for that in most cases. Never- 
theless, in order to make any conclusions about the rate 
of microvariability, the influence of this and other pos- 
sible effects has to be fully quantified with the proper 
efficiency calculation. We present such a calculation in 
Section EQ1 

The approach taken by dL98 in the Hipparcos analysis 
is not possible for OGLE-II data. Hipparcos light curves 
consist of brief observing sequences lasting typically less 
than 1 day but densely sampled with measurements often 
separated by about 20 minutes. This is why long term 
variations of the main Mira cycle are not an issue and at 
the same time many events can be resolved. The time 
sampling of the OGLE-II light curves is generally not 
sufficient to resolve any rapid variability on time-scales 
shorter than 1 day. Typically the density of measure- 
ments well within the observing season corresponds to 
one point every 1 to 3 days and it tapers off near both 
ends of the season. Some variations of temporal sampling 
arise due to changing observing conditions. Still, there 
are 1046 pairs and 210 triples of observations spanning 
less than 24 hours. Those better sampled regions sig- 
nificantly improve the sensitivity to the shortest events 
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Fig. 2. — Examples of OGLE-II light curves for Mira variables in the Galactic Bulge with model fits given by Equation Q (left). Also 
shown are photometric residuals with respect to the model (right). All units are / magnitude versus time in days. Note that the amplitudes 
of all events reported by dL98 exceed the largest residuals by a factor of ~5. 



considered in our analysis. At the basis of our approach 
is exceptionally good error distribution in the OGLE-II 
DIA photometry combined with the ability to fit longer 
term behavior of the Mira flux. Even if microvariability 
events are infrequent and very short, each object with 
rapid variations discovered by dL98 spends some (small) 
fraction of time in the events, when the object flux is 
far from the normal smooth trend of the main pulsation 
cycle. Eventually, with persistent monitoring, some mea- 
surements will be affected by microvariability regardless 
of the time sampling. 

We define an "event" as one or more consecutive mea- 
surements deviating (consistently) up or down by more 
than 0.2 mag with respect to the best fit model of Equa- 
tion^ Experiments requiring two consecutive deviations 
(a two-point filter) and 0.15 mag limit for the flux depar- 
ture generally confirm our conclusions, albeit at a factor 



of 2 to 10 lower sensitivity. In both searches we detect 
no events. 

We checked the impact of the last selection cut from 
Section 12.31 which removes points with unusually large 
error bars. Without this condition imposed, only three 
single point "events" in the entire data set depart from 
the model by: 0.21, —0.25, 0.23 mag. Two of them orig- 
inate from the same frame, and all have error bars much 
larger than the rest of the corresponding light curve, par- 
ticularly the neighboring points. They are much more 
likely manifestations of a problem with the measurement 
rather than detections of microvariability. It is safe to 
conclude that we have not detected any rapid variability 
events. 

3.2. Probability density for event rate 

Counting events from Section l3~T1 is a Poisson process 
with the expectation value = ^ s vTe s — v£, where v 
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is the yearly rate of actual physical events per object of 
the ensemble, T is the duration of the survey and e a is the 
efficiency of detecting an event that occurred randomly 
during the survey time in star s. We substitute £ = 
T ^2 S e s to simplify notation. In other words fi is the total 
expected number of events in the entire sample of 485 
stars observed for T = 3 years assuming that every star 
undergoes v events per year on average. Equivalently, if 
only a fraction of the stars experienced events, the rate v 
would be proportionately higher in those objects. Using 
Bayes' theorem we can derive the probability distribution 
for the rate given the number of observed events n: 



P{y\n) 



P(n\v)P(v) 



(2) 



where P(n\v) = tf 1 exp(— fJ,)/n\ is the Poisson distribu- 
tion and P{n) is the normalization factor obtained from 
condition J °° P(y\n)du — 1. For lack of other knowl- 
edge we also assume a uniform prior probability for the 
rate P{v) = 1. This results in the following probability 
distribution for v; 



P(y\n) 



£ 



n+l 



-v n exp(— 1/£). 



(3) 



In the case of a non-detection P(v\n = 0) = £exp(— v£) 
which can be integrated for a confidence level a to obtain 
the upper limit for the rate: 



ln(l - a) 
£ ' 



(4) 



In the following discussions we will use i?g 5 , the 95% 
confidence upper limit for the event rate. 

3.3. Simulating efficiency 

The term £ is the crucial piece of information and can 
be evaluated with the use of a Monte Carlo simulation. 
We consider microvariability events with two kinds of 
flux profiles: triangular and instantaneous rise with ex- 
ponential decay (Sections 14. II and 14.21 respectively) . The 
time-scale (r) and the amplitude of a simulated event 
(Am) are drawn from their assumed distributions. A 
random epoch within the three-year survey baseline is 
also generated and the flux variation is injected at that 
epoch. Finally, we attempt to detect the event using the 
same procedure employed with the real observations in 
Section 13.11 This includes re-fitting the empirical curve 
describing normal smooth behavior. For each of the 485 
stars in the ensemble, we generate 1000 trial events and 
count all successful detections. 

4. RESULTS 

4.1. Derived rate of Hipparcos events 

The distributions of time-scales and amplitudes for mi- 
crovariability events reported by dL98 are shown in Fig- 
ure For the purpose of constraining the rates of these 
types of event, the simulated event time-scales and ampli- 
tudes are drawn from the actual observed distributions. 
There were approximately as many upward and down- 
ward deviations reported in the Hipparcos data. We ver- 
ified that the two distributions appear uncorrelated. In 
all our simulations we assume that the event amplitude 




Fig. 3. — Distribution of time scales (left) and amplitudes (right) 
of events reported by de Laverny et al. (1998). 



is independent of the time-scale and the sign of the de- 
viation is random. The shape of the artificial impulse is 
a triangle with the base covering the time interval ±r. 
In Section T4.2I we will see that the results are not very 
sensitive to that choice. 

Figure 0] shows the probability distributions for the 
event rate v given n = 0, 1, 2, and 3 detections. Given 
that we have detected n = events in the OGLE-II data, 
the upper limit for the rate (95% confidence) is -R95 = 
0.038 yr _1 per star. We conclude, with 95% probability, 
that dL98 type events are less frequent than 1 every 26 
years in any given object of our OGLE /-band sample. 

The study of dL98 does not give the value for the ef- 
ficiency of the Hipparcos light curves in detecting Mira 
microvariability. In order to derive the event rate im- 
plied by results of dL98 we need to know the detection 
efficiency. We can estimate it using public Hipparcos 
data from the Epoch Photometry Annex (ESA 1997). 
For any given object the Hipparcos observing sequences 
cover only a small fraction of the total duration of the 
experiment. Observing sequences are typically separated 
by gaps lasting about two weeks. Any events occurring 
in those gaps cannot be detected; the instrument is sim- 
ply not observing the object. We extracted light curves 
for 226 Hipparcos stars classified as Miras in GCVS. This 
is fewer than 239 light curves analyzed by dL98, but in- 
cludes objects contributing 50 out of 51 reported events 2 . 
We adopt 3.37 years as the total duration of the Hippar- 
cos mission, taking the difference between the first and 
the last epoch in all 226 light curves. For each light curve 
we subtract from the total duration the sum of all gaps 
lasting 10 days or longer and divide the result by the 
total duration. On average there were 104 photometric 
points available per star. We do not remove any flagged 
points. This and the implicit assumption of 100% effi- 
ciency within observing sequences will only increase the 
final value of the efficiency and lower the rate. 

The distribution of the resulting efficiencies for Hippar- 
cos data is shown in Figure[S]and compared to the results 
of our Monte Carlo experiments with OGLE-II data. For 
a few Hipparcos stars events can be detected with rela- 
tively high efficiency reaching 20-40%. However, for the 
majority of the Hipparcos light curves the efficiency is 
similar to that for OGLE-II data. The additional Monte 
Carlo run on OGLE-II data with the distribution of time- 
scales from dL98 truncated at r > 2 days demonstrates 
that our results do not depend dramatically on the pres- 
ence of the long duration tail. 

If v is the true yearly rate of events in units of yr _1 

2 There seems to be a problem with identification of the 51st 
event in variable CE Lyr. The on-line Hipparcos catalog entry for 
the star closest to the GCVS position states that the star may be 
wrongly identified with CE Lyr. The star appears constant. In any 
case, these minor differences are unimportant. 
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Fig. 4. — Probability distributions for the rate of events of the de 
Laverny et al. (1998) type. Lines correspond to probabilities given 
that there were between n = (farthest left) to n = 3 (farthest 
right) detections. The location of the vertical line indicates the 
95% confidence upper limit for the rate given a non-detection. 
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star -1 , the expected number of events over the total 
duration of the experiment (detected or not) is: vNT, 
where N is the total number of monitored objects and T 
is the duration of the experiment in years. On the other 
hand, the total number of events can be recovered from 
the number of detected events by counting each detec- 
tion with the weight proportional to the inverse of the 
efficiency: n s/^s, where n s is the number of detected 
events in star s and e s is the efficiency of detecting a ran- 
dom event in star s. The number of detections per star 
in dL98 study is between and 3. By equating the two 
above expressions one obtains an estimator of the event 
rate: 

s 

We evaluate the rate in Equation[S]using our efficiency 
estimates for a sample of N — 226 Hipparcos light curves 
and accepting all detections reported by dL98 (except the 
one in CE Lyr for which we could not find the data). As 
mentioned earlier, the efficiency was likely over-estimated 
so the rate is likely under-estimated. The resulting rate is 
1.63 yr _1 star" 1 . Therefore, the Hipparcos results imply 
that microvariability events in Miras are quite common 
in the spectral region covered by the Hp band. This is 
due to the fact that, assuming the events are not corre- 
lated with the scheduling of the Hipparcos observing se- 
quences, for every detected event there were ~25 events 
that could not be detected due to a gap in coverage. The 
effective efficiency is only ^4% and comparable to that 
of the OGLE-II observations. Taken at face value, the 
OGLE-II data appear inconsistent with high rates im- 
plied by the number of Hipparcos detections. Given a 
true rate of 1 event per year per object, the probability 
that an ensemble like ours would yield detections is less 
than 6 x 10 -35 . However, there are several scenarios that 
can easily accommodate both our non-detection in the I 
filter and the positive Hipparcos result in the Hp band. 
We discuss this issue further in Section |SJ 



4.2. Constraining flashes on Miras 

Several possible mechanisms for rapid light variations 
in Mira variables have been discussed by Willson and 



Struck (2002). One such scenario invokes a sudden re- 
lease of energy leading to " hot flashes" . We can put use- 
ful constraints on such events by simply simulating the 
efficiency of detecting a different kind of flux impulse in 
our data. This also provides a cross-check on the results 
from Section 14.11 and in particular shows how sensitive 
the derived limits on frequency are with respect to the 
pulse shape. In this section we adopt a flare with in- 
stantaneous rise by Am and exponential decay time r. 
The average efficiency (per star) of detecting such flares 
is shown in Figure The dependence on the flare am- 
plitude between 0.3 and 1.1 mag is plotted for a number 
of time scales r ranging from 0.1 to 1.6 days. Note that 
the figure shows that flashes with amplitudes of a few 
tenths of a magnitude lasting only a couple hours can be 
detected with the efficiency of a few percent in any of the 
stars in our ensemble. 

Before deriving the upper limits on rates, we intro- 
duce the size-frequency distribution into this calculation. 
The published data offers very little information on the 
properties of suspected flares, but one generally expects 
strong events to occur less frequently and weak events 
to dominate. For the purpose of illustrating the range of 
possible outcomes, we assume a power-law distribution 
of / = AF/F the amplitude of the simulated flashes rel- 
ative to the flux F of the parent star: 

N(f)df cx ndf (6) 

At fixed time scale r and parent star flux F this size dis- 
tribution is equivalent to assuming a power-law distribu- 
tion of flare energies. The distribution is truncated be- 
low the relative flare amplitude of 20% and above 100%, 
which corresponds to amplitudes of about 0.2 and 0.75 
mag respectively. 

Taking the amplitude distribution of dL98 and ignoring 
the sign of the variation, 7 between 2.0 and 2.5 gives a 
good approximation within the limited statistics of the 
data set. We evaluate the efficiency at several fixed time 
scales ranging from 0.1 to 2.0 days for a power-law size 
frequency distribution with index 7 =1, 2 and 3. The 
resulting rates are given in Figure The frequency of 
1 per 2.5 years in any given star is excluded with 95% 
probability for all types of considered events down to 
durations of 2.4 hours. The limits quickly tighten for 
longer flare durations. For 7 = 2 we can exclude 0.5 day 
long events being more frequent than about once every 
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Fig. 6. — Simulated efficiency of detecting flashes with instanta- 
neous rise and exponential decay. The dependence on amplitude 
Am is given for several fixed time scales t. 




Fig. 7. — Upper limits (95% confidence) on the rate of flashes 
on Mira variables given a non-detection in our ensemble. The de- 
pendence on time scale is given for three power law distributions 
of fr actio nal flux amplitudes with the index 7 = 1,2, and 3 (Sec- 
tion 



17 years. 

5. DISCUSSION 

Using a large sample of objects with high quality light 
curves we place stringent limits on the frequency of rapid 
variability events in the Galactic Bulge Mira variables in 
the / band. To the best of our knowledge this is the first 
published limit for the rate of rapid variability in Miras. 
The study is most sensitive to events lasting about 1 day 
or longer, however it retains a few percent efficiency (per 
star) for detecting unresolved microvariability on time 
scales as short as 2 hours. With 95% probability we ex- 
cluded events of the type found by dL98 occurring more 
frequently than 1 per 26 years in the average object of our 
sample (0.038 yr _1 star -1 ). The results are not sensitive 
to the details of the temporal profile of the flares. The 
most significant dependence of the detection efficiency 
in our ensemble is on the event time scale because of 
the 1-3 day light curve sampling. The limits on the fre- 
quency of flashes with a sudden flux rise and exponential 
fall-off with power-law distribution of energies (at fixed 
parent star luminosity) have been presented. For flare 
durations ranging from 0.1 to 1 day the corresponding 
limiting rates are between once every few years and once 
every few tens of years. 

In our calculations we essentially extrapolate the prop- 
erties of the rapid variability events observed by Hippar- 
cos from Hp band to the / band. Under that assump- 
tion, the OGLE-II data are inconsistent with high rates 
implied by the number of Hipparcos detections. A possi- 
ble explanation is that the microvariability is much less 
prominent or the events are shorter in the I band. In that 
case we are still learning something about microvariabil- 



ity in Mira stars. The typical accuracy of the OGLE-II 
photometry is much better than 0.2 mag (Figure EJ) and 
we could not find any indication of rapid variations in 
any of the 485 studied objects. All light curves have 
very smooth appearance. 

There are at least three possible explanations for the 
difference between the two results: 

1. Mira variables in the Galactic Bulge population dif- 
fer from their counterparts in the solar neighbor- 
hood. In the scenario of a planet interacting with 
the Mira atmosphere, e.g., the Bulge population 
might not have any planets. It is plausible that the 
Galactic Bulge Miras have lower metallicity than a 
typical Hipparcos star. Planets are more prevalent 
around metal-rich stars (e.g. Santos, Israelian, & 
Mayor 2004) . It has been suggested that dust may 
play an important role in microvariability of Mi- 
ras (Stencel et al. 2003). The dust content of the 
stellar envelope could be related to metallicity. 

2. Microvariability in Miras is occurring in the part 
of the spectrum that is not covered by Cousins I 
band but is contributing to the blue part of the Hp 
band. 

3. Microvariability in Miras is less common than im- 
plied by results of dL98. The only sure way to ver- 
ify this is to collect more data and carefully analyze 
it. Numerous groups are setting up fast cadence 
photometric monitoring experiments with the goal 
of detecting real time transients (Vestrand et al. 

2002) or finding transiting Jovian planets (Home 

2003) . Teets & Henson (2003) are conducting a 
special purpose monitoring campaign of ~5 sources 
previously reported to display rapid variability us- 
ing B, V, R, I filters. The current data set covers 
many nights spanning the time interval of over 2 
years with up to 5 photometric points per hour for 
up to 8 hours a night. So far only a single event 
in RR Boo has been detected. With the available 
phase coverage this is far below expectations based 
on the Hipparcos data. 

The Hp filter of the Hipparcos mission is broader than 
V and still sensitive in the wavelength range 650-900 
nm, where the spectral energy distribution of a Mira 
variable rises steeply. We can calculate the amount of 
energy emitted by a Mira star in the Hp band know- 
ing that the ratio of luminosities in Hp and V is about 
2.5. This is easily estimated assuming that the spectrum 
of a Mira star blue-ward of 1000 nm is a black body 
with effective temperature T e // = 2500 K and by fold- 
ing it with filter transmission curves (Johnson 1965, van 
Leeuwen et al. 1997). Then using the absolute magni- 
tudes My = —0.3 of an M5 giant and the U-band Solar 
luminosity Ly.Q — 4.64 x 10 32 erg s _1 (Binney & Mer- 
rifield 1998) we obtain L Hp = 1.3 x 10 35 erg s _1 for a 
Mira star. We adopt My iQ = 4.83. Mira variables vary 
considerably less in the Hp band compared to V with 
amplitudes rarely exceeding 4 mag (Whitelock, Marang 
& Feast 2000). Even during minimum light and taking 
an extreme value for the amplitude (4 mag) the total 
energy of a flare with amplitude Am and duration r is: 
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the same size. The observational evidence is accumulat- 
ing and it should not be long before we can formulate a 
(7) more definitive answer to the exciting question of rapid 
high amplitude variability on Miras. 



where C is the numerical factor depending on the profile 
of the light variation. For a top hat function C = 1. We 
can see that a flare of a third of a magnitude lasting 24 
hours will require 10 38 ergs of energy — an order of magni- 
tude more than the energy required for the same contrast 
in the V band calculated by Willson & Struck (2002). If 
in fact the events detected by Hipparcos are confirmed, 
one will be faced with finding a large energy source that 
is triggered approximately as often as an energy sink of 
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